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ABSTRACT 

Weld  properties  in  Friction  Stir  Welding  (FSW)depend  on  the  peak  temperature  attained  during  the  process  and 
cooling  rate.  Temperature  cycles  in  this  joining  process  decides  the  microstructure  not  only  in  the  weldment  but  also  in 
the  heat  affected  zones  of  the  parental  metal.  The  complex  flow  of  plasticized  metal  makes  it  difficult  to  measure 
generated  heat  directly  which  results  the  metal  mixing  and  flow  in  the  stir  zone  in  FSW  process.  This  creates  an 
importance  on  the  analytical,  numerical  and  computational  modelling  of  friction  stir  welding.  In  this  review  on  friction 
stir  welding,  special  focus  is  given  on  various  analytical  models  developed  to  estimate  peak  temperature  during  the 
process.  The  different  approaches  that  have  been  used  to  investigate  heat  generation  are  presented  and  their  drawbacks 
are  discussed.  Assumptions  made  in  the  modelling  and  their  corresponding  effects  in  the  results  are  also  addressed. 
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INTRODUCTION 

Friction  Stir  Welding  (FSW)  is  one  of  the  metal  joining  processes,  performed  in  plasticised  solid  state 
parental  metal  with  the  help  of  frictional  heat.  In  a  metal  joining,  the  weld  properties  are  influenced  by  the  peak 
temperature  during  the  process  and  rate  of  cooling  after  the  joining  process  [1],  FSW  produces  comparatively 
better  weld  properties  as  it  is  performed  in  much  lower  temperature  than  other  conventional  joining  processes  [2]. 
This  solid  state  metal  joining  process  is  performed  by  forcing  a  non-consumable  rotating  tool  to  penetrate  and 
move  through  the  entire  joint  of  the  parts  to  be  welded  (Figure.  1).  All  the  welding  parameters  in  FSW  are  selected 
to  soften  the  work  piece  material  to  facilitate  mechanical  deformation.  The  frictional  heat  produced  in  the 
tool/matrix  interface  plasticises  and  joins  the  parental  metal  in  the  line  of  joint  and  the  tool  [3]. 

FSW  operations  can  be  divided  into  several  phases  based  on  the  relative  motion  absorbed  between  work 
piece  and  tool  [5].  Initially  tool  pin  is  rotated  and  allowed  to  pierce  vertically  into  the  work  piece  called  plunging 
operation  and  the  material  is  too  cool  for  plastic  deformation.  In  this  stage,  rubbing  action  created  by  the  pin 
creates  chips  like  other  machining  process  and  temperature  is  increased.  When  the  tool  shoulder  touches  the  weld 
surface  the  vertical  movement  is  stopped  and  tool  is  allowed  to  rotate  for  a  while  called  dwell  period.  The  entire 
tool  shoulder  and  pin  surface  contribute  to  the  frictional  heating  and  the  work  piece  reaches  a  critical  temperature 
for  plastic  flow.  After  these  two,  welding  operation  is  carried  out  by  the  horizontal  movement  of  the  welding  tool 
along  the  joint  line.  Most  of  the  welding  defects  like  porosity  and  void  formation  in  the  conventional  joining 
processes  are  eradicated  in  FSW  as  the  shoulder  confines  the  underlying  material.  During  the  plunge  action,  the 


www.tjprc.org 


SCOPUS  Indexed  Journal 


editor@tjprc.  org 


Original  Article 


180 


J.  Stephen  Leon  &  V.  Jayakumar 


conical  shoulder  acts  like  an  escape  volume  of  the  material  displaced  by  the  tool  pin  [6]. 


Figure  1:  Different  Phases  in  Friction  Stir  Welding  Process  [4] 

Almost  all  the  properties  of  the  final  weld  depend  on  generating  heat  and  heat  flow  through  it  during  the  joining 
process  [7].  So  understanding  the  temperature  profile  during  the  process  is  important  to  optimise  the  weld  tool  geometry 
and  other  weld  parameters  in  order  to  obtain  improved  weld  microstructure  and  other  mechanical  properties  associated 
with  it. 


Tool  Shoulder 


Figure  2:  Different  Zones  in  the  FSW  Process  [4] 

This  thermo  mechanical  stir  metal  joining  process  produces  two  different  zones  namely  stir  zone  and  heat  affected 
zone  (Figure. 2).  Intense  deformation  produced  by  the  FSW  tool  rotation  makes  all  conventional  temperature  measurement 
techniques  like  embedding  thermocouples  in  the  stir  zone  impossible  [8].  This  necessitates  an  analytical  thermal  modelling 
to  predict  generated  heat  to  define  thermal  profile  during  friction  stir  welding  process.  Therefore  the  maximum  temperature 
is  estimated  indirectly  by  analysing  the  change  in  the  microstructure  like  the  dissolution  of  larger  precipitates  [9]  in  the  stir 
zone  of  the  weld  or  by  measuring  the  temperature  in  the  region  adjacent  to  the  stir  zone  with  the  help  of  embedded 
thermocouples  and  predicting  the  peak  temperature  in  the  stir  zone  by  analytical  methods  as  the  temperature  is  a  function 
of  distance  [10].  These  indirect  estimations  couldn’t  produce  the  promising  result  as  sometimes  the  predicted  temperatures 
are  more  than  melting  point  [11]  which  is  not  true  in  FSW  as  it  is  a  solid  state  joining  process.  This  is  due  to  the  fact  that 
this  increase  in  the  dislocation  density  and  an  increase  in  the  amount  of  the  grain  surface  and  grain  edge  per  unit  volume  is 
not  only  because  of  a  temperature  rise,  but  also  due  to  the  plastic  deformation  produced  by  the  tool  during  stirring  action 
[12]. 

In  FSW,  heat  is  generated  by  the  friction  developed  in  the  tool  workpiece  interface  and  also  due  to  the  plastic 
deformation  that  occurs  away  from  the  tool  workpiece  interface.  Local  heat  generation  produced  by  the  deformation  and 
grain  boundary  sliding  contributes  only  4.4%  of  the  total  heat  produced  [13]  and  the  remaining  heat  generation  is  because 
of  frictional  heat.  So  in  this  paper,  special  emphasis  has  been  given  on  the  analytical  estimation  of  heat  generation  due  to 
the  friction  in  the  tool/matrix  interface. 
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Thermal  modelling  is  used  to  analyse  temperature  distribution  during  FSW.  All  thermal  models  in  FSW  are 
interested  to  solve  a  general  heat  conduction  equation  (1),  by  opting  different  initial  condition  and  appropriate  boundary 
conditions  [14]. 


pCpT  +  pfiCpTi  -  ( kTi)  i  + 


(1) 


Here  Qgen/V  denotes  volumetric  heat  source.  Unlike  other  joining  process,  in  FSW  heat  is  generated  not  only  by 
means  of  external  sources,  but  also  generated  by  friction  and  plastic  dissipation  in  the  tool/matrix  interface  [15]. 
Tang  et  al.[16]  analysed  the  tool  shoulder  and  the  tool  pin’s  contribution  to  the  heat  generation  and  concluded  that  as  the 
shoulder  has  contact  area  and  contact  pressure  comparatively  much  higher  than  the  pin,  shoulder  dominates  the  heat 
generation.  Assuming  only  sliding  condition  and  heat  generation  by  the  shoulder  in  the  interface,  a  simplified  expression 
for  heat  generation  (eqn2)  is  derived  by  implementing  numerical  model  based  on  position  dependent  surface  flux  for  flat 
shoulder  FSW  tool  neglecting  the  heat  generated  by  the  pin[17]. 


Rtotal  — 


3 Qtotalr 
2nRihoulder 


(2) 


This  is  a  commonly  used  expression  to  calculate  heat  generation.  In  this  equation  Q  totai  is  an  input  parameter. 
Unlike  other  metal  joining  process,  Q  totai  is  not  controlled  by  external  sources.  But  it  itself  depends  on  internal  function 
variables.  So  this  equation  has  limitations  and  cannot  be  used  to  predict  all  FSW  welding  parameters.  For  example,  Q  totai 
itself  is  a  function  of  rotational  speed  [18].  In  most  of  the  thermal  modelling,  Q  total  is  measured  using  dynamometer  and  is 
taken  as  a  constant  to  find  q  total[19]. 


Thermal  model  is  used  to  analyse  temperature  distribution.  In  FSW,  as  generated  heat  is  a  function  of  frictional 
contact  conditions  and  material  flow  it  cannot  be  a  constant  value  [20].  Heat  generation  also  depends  on  tool  geometry  and 
other  weld  parameters  [18].  So  pure  thermal  model  has  limitations  and  is  replaced  by  thermo  mechanical  or  CFD  model 
which  couples  heat  generation  and  material  flow  by  plastic  dissipation  in  the  shear  layer  and  frictional  contact  at  the  tool 
and  work  piece  [21]. 


THERMAL  MODELING  BASED  ON  THE  SOLIDUS  TEMPERATURE 


It  was  observed  from  the  various  experimental  analysis  [22-25]  that  peak  temperature  during  FSW  lies  between 
60  to  90%  of  the  solidus  temperature  of  parental  metal.  Arbegast  and  Hartley  [26]  derived  an  expression  (eqn3)  which 
relates  the  peak  temperature  developed  during  the  process  and  solidus  temperature  of  the  parental  metal. 

—  =  (  0,2  )  (3) 

Tm  Vv  X  10 V  v  ; 

where  a  is  in  the  range  of  0.04  to  0.06,  k  is  the  thermal  conductivity,  co  and  V  are  tool  rotation  and  traverse  speed. 

During  the  joining  process  if  the  generated  temperature  approaches  solidus  temperature  of  the  material  to  be 
joined,  the  material  in  the  tool  metal  interface  will  act  like  a  fluid,  which  reduces  friction  in  the  tool  metal  contact  surface 
tremendously  and  results  ina  negligible  heat  generation  and  temperature  brought  down.  This  self-stabilizing  effect  always 
maintains  the  temperature  below  its  solidus  temperature  [26].  In  an  approach  developed  by  Tutum  et  al.[27]  peak 
temperature  in  the  tool  metal  interface  is  taken  as  a  constraint  to  optimize  other  weld  parameters.  This  approach  helps  to 
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model  temperature  profile  without  prior  knowledge  about  heat  input.  A  thermal  model  developed  based  on  this 
assumptions  deviates  from  actual  conditions  as  the  peak  temperature  is  not  constant  throughout  the  process. 
Varying  temperature  during  the  process  produces  a  non-uniform  thermo  mechanical  conditions  which  are  not  considered  in 
this  thermal  modelling.  A  new  thermo  pseudo  mechanical  formulation  can  overcome  these  limitations  which  consider  the 
varying  conditions  in  the  prediction  of  heat  generation  [28]. 

THERMO  PSEUDOMECHANICAL  MODELLING 

In  this  approach,  the  total  amount  of  heat  generation  is  modelled  from  frictional  heat  and  heat  developed  due  to 
the  plastic  dissipation  at  the  tool  metal  interface  [17-19].  Colomb’s  law  of  friction  is  used  to  determine  the  shear  stress 
between  the  contact  surfaces  (eqn4).  A  uniform  thickness  shear  layer  was  prescribed  analytically.  The  material  flow  in  the 
shear  layer  produces  convective  heat  generation  in  the  weldment,  around  the  probe. 

r  =  pP  —  pa  (4) 

Hasimoto  et  al.[29]  observed  that  the  peak  temperature  during  FSW  is  a  strong  function  of  tool  traverse  and 
rotation  speed.  Based  on  the  influence  of  contact  pressure  and  tool  rotation  rate,  Friggard  et  al[30]  developed  a  model 
(eqn5)to  calculate  heat  input  per  unit  area  as 

q0  =  ^ n2pPtoR 3  (5) 

Where  p  is  friction  coefficient,  P  is  contact  pressure,  at  angular  velocity  of  tool  and  R  is  radius  of  tool  shoulder. 

Heat  generation  due  to  the  friction  is  the  product  of  the  tangential  speed  of  the  tool  with  respect  to  the  workpiece  at  the  tool 
work  piece  interface  and  the  frictional  force  and  given  by  (eqn6)[31], 

dQ  =  (cor  —  Vsin  6)ppdA  (6) 

Where  p  is  local  pressure  and  d  A  is  element  area. 

Normally  heat  is  generated  in  the  surface  very  close  to  the  tool  and  has  a  complex  geometry  based  on  the  tool 
geometry.  But  for  analytical  estimation,  tool  design  is  assumed  to  be  conical  shoulder  with  a  cone  angle  a  or  for  flat 
shoulder  surface  in  which  a  =  0.  VsinB  can  be  neglected  as  the  effect  of  tool  rotation  speed  (cor)  is  much  higher  than  the 
effect  of  tool  traverse  speed  (VsinB)  [32&33]  so  heat  generation  equation  can  be  simplified  as  (en7) 

dQ  =  codM  —  cordF  —  a>rTcontactdA  (7) 

In  equation  (7),  for  dA  is  achieved  considering  circular  surface  perpendicular  to  the  rotation  axis  for  probe  tip, 
cylindrical  surface  parallel  to  the  rotation  axis  for  probe  side  and  conical  surface  tilted  with  respect  to  the  rotation  axis 
(Figure. 3). 
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Figure  3:  Heat  Generation  Contributions  in  Analytical  Estimation  [36] 


Figure  4:  A  Schematic  Drawing  of  Surface  Orientations  and  Infinitesimal 
Segment  Areas  for  Probe  Tip,  Probe  Side  and  Conical  Tool  Shoulder  [36] 


Total  heat  generated  is  an  addition  of  heat  generated  under  the  tool  shoulder  (Qi),  tool  probe  aside  (Q2)  and  probe 
tip  (Q3)  (Figure  4).  Heat  generation  equations  for  QhQ2  and  Q3  are  derived  from  general  heat  generation  equation[34]. 
Adding  all  final  equation  derived  is  (eqn8). 


Q total  3  tTT contact^  . 


shoulder  ~F  ^^-probe^probe) 


(8) 


In  this  equation  shear  stress  value  is  kept  constant  in  the  contact  interface.  But  in  reality,  the  contact  shear  stress 
varies  depending  on  contact  conditions  [35].  To  define  contact  yield  stress,  three  different  fully  sticking,  fully  sliding  and 
partial  sticking/sliding  contact  conditions  have  to  be  analysed  in  the  tool/matrix  interface  [37]. 


Sticking  Condition 

The  work  piece  is  divided  as  surface  segment  and  lower  surface  segment.  If  the  shear  yield  stress  due  to  the 
frictional  contact  between  the  tool  and  the  work  piece  is  more  than  the  yield  stress  of  the  underlying  segment,  it  results  in 
full  sticking  condition.  In  this  condition,  the  matrix  segment  accelerates  till  internal  matrix  shear  stress  and  contact  shear 
stress  become  equal  and  finally  reaches  tool  velocity.  Schmidt  et  al  [36]  defined  a  contact  state  variable  (5)  (eqn9)  as 

g  —  Vmatrix  _  ^  _  Y  /m 

Ytool  ^  tool 


where  y=  slip  rate  and  V  tooi  =  rco. 

During  sticking  condition  (5=1  (or)  V  tooi  =  V  matrix).  Shear  stress  occurs  in  the  layer  very  close  to  the  tool/matrix 
interface.  In  this  case,  there  exists  a  change  in  velocity  from  the  layer  sticking  with  the  tool  till  the  stationary  layer  in  the 
matrix.  So  the  deformed  layer  starting  from  the  tool  interface  should  be  treated  as  a  shear  surface.  So  shear  stress  x  contact  = 
x  yield  =  o  yieid  /V3.  The  extent  of  slip  can  also  be  estimated  by  experimentally  measured  values  using  the  expression  (eqnlO) 
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given  by  Deng  et  al  [38]. 

S  =  0.2  +  0.8  [l  -  exp  (-S0  ^)]  (10) 

Where  at  r  represents  the  torque  at  the  point  of  calculation  and  co0Rs  is  the  maximum  torque  that  can  be  produced 
by  the  tool. 

Sliding  Condition 

Sliding  condition  occurs  when  friction  shear  stress  is  lesser  than  internal  matrix  yield  stress.  In  this  condition 
matrix  segment  slightly  elastically  deforms  and  shear  stress  is  equal  to  the  dynamic  contact  shear  stress.  During  fully 
sliding  condition  (8  =  0  (or)  V  matrix  =  0).  Shear  stress  occurs  in  the  contact  interface.  Here  r  contact  =x  friction  =  qP  =  qo.  Here 
P  and  a  are  contact  pressure.  In  a  thermal  model  developed  by  Hilgert  et  al. [39]  pure  sliding  condition  is  adopted  and  heat 
generation  equation  is  simplified  as,  q  total  =  co  tool  r  x  yield 

Partial  Sticking/Sliding  Condition 

In  this  contact  shear  stress  is  equal  to  the  internal  matrix  yield  stress  due  to  a  quasi-stationary  plastic  deformation 
rate.  The  matrix  segment  accelerates  and  stabilises  in  a  velocity  lesser  than  the  tool  surface  velocity  n  a  new  approach, 
Schmidt  et  al.[19]  derived  an  equation  (eqnll)  considering  frictional  and  plastic  dissipative  heating  under  varying  contact 
yield  stress. 

Rtotal  —  YT friction  "P  (<Wooi^  ~  Y^yield  (H) 

where  q  totai  is  the  total  heat  flow,]/  is  the  shear  rate,  TfriCtion  is  the  shear  stress  due  to  friction,  x  yie]d  is  shear  yield 
stress,  co  tooi  is  the  angular  velocity  of  the  tool  and  r  is  the  distance  from  the  heat  source  to  the  centre. 

An  energy  dependent  expression  developed  by  Hamilton  et  al.  [40]  reveals  thatslip 

8E  =  exp(-^-)  (12) 

'  v^Umax' 

when  (E^eff  =  0  then  5  E  =  1  and  sticking  friction  dominates.  When  (E,)  elT  increases,  metal  in  the  tool/metal 
interface  softens  and  slip  occurs.  So  maximum  possible  case  can  be  (E|)rff  =  (E,)[Ilax  and  in  this  case  5  E  =  0.37.  From  this  it 
can  be  concluded  that  slip  rate  can  lies  betweenl  to  0.37. 

EXPERIMENTALLY  FOUND  TORQUE  BASED  THERMO 
PSEUDOMECHANICAL  MODELING 

In  the  modelling  of  temperature  profiles  Rosenthal  solution  is  used.  The  input  power  during  friction  stir  welding 
process  has  been  used  as  input  in  the  various  models  [41,  42]  developed  based  on  Rosenthal  equation.  Khandkar  et  al.  [43] 
developed  heat  input  models  based  on  torque  input  as  the  heat  input  is  equal  to  the  rotational  power.  This  experimentally 
arrived  torque  value  replaces  the  heat  input  to  find  yield  stress  on  the  contact  surface  in  the  elevated  temperature. 
This  analytical  solution  created  a  link  between  plunge  force  and  heat  generation  as  the  plunge  force  is  used  to  calculate 
torque  input. 

T total  =  2qF(J+|h)  (13) 

where  F  is  applied  force,  r0  is  shoulder  radius,  r,  is  radius  of  pin  and  h  is  pin  height.  The  obtained  analytical  results 
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compared  with  experimental  obtained  by  Frigaard  et  al  [44  &4  5]  (Figure  5)  to  analyse  the  relationship  between  plunge 
force  and  heat  generation. 


Figure  5:  Comparison  of  Plunge  Force  Experienced  [44] 

During  plunge  period,  plunge  force  and  torque  increases  steadily.  During  dwell,  in  analytical  solution  it  is 
predicted  to  be  sliding  and  uniform  pressure  condition.  But  from  the  experimental  values  it  is  observed  that  plunge  force 
and  torque  reduces.  During  welding  period,  little  increase  in  the  plunge  force,  but  torque  remains  the  same. 
Thus  an  increase  in  the  plunge  force  does  not  affect  torque.  It  is  understood  that  during  welding,  a  reaction  force  is 
developed  which  opposes  plunge  force  and  tries  to  move  the  tool  upwards.  To  overcome  the  reaction  force,  the  machine 
gives  extra  force.  Thus  at  this  stage,  plunge  force  cannot  be  used  to  access  the  heat  generation. 

Power  input  given  to  the  machine  is  influenced  by  tool  rotation  speed  (to)  and  tool  transverse  speed  (v)  [46]  and  is 

given  by 

P  tool  =  M  a>  +  Fx  V  (14) 

Where  F  x  is  forward  force,  M  denotes  torque.  The  relationship  between  torque,  tool  rotation  speed  and  tool 
traverse  speed  is  analysed  [47  &48]  and  it  was  concluded  that  torque  could  be  reduced  by  the  increase  in  tool  rotation  and 
reduce  in  tool  traverse  speed.  The  influence  of  tool  rotation  speed  is  higher  than  tool  traverse  speed.  Based  on  this  an 
equation  was  developed  by  L  .  Cui  et  al.  [15]  which  relates  the  input  power  and  tool  rotation  speed 

M  =  M0  +  Mf  e"n“  (15) 

where  M0  is  the  minimum  torque  required,  Mf  is  the  pre  exponential  parameter  and  n  denotes  the  decay  parameter. 
Upadhyay  and  Reynold’s  [41]  study  on  Peak  temperature  (T)  during  FSW  with  respect  to  tool  rotation  (<n)  reveals  that  T  is 
a  function  of  <n  and  they  expressed  peak  temperature,  based  on  the  graph  drawn  between  T  and  co,  as 

T  =  T0  +  T  f  (l-e"q“)  (16) 

THERMO  PSEUDOMECHANICAL  MODELING  CONSIDERING  SURFACE  FLUX 

Yield  stress  in  the  shear  layer  is  independent  of  pressure.  But  it  is  highly  influenced  by  temperature. 
In  the  previous  models  used  to  find  heat  generation,  in  order  to  apply  constant  shear  yield  stress,  an  isothermal  is  assumed 
in  the  tool/matrix  interface.  But  in  real  practice,  heat  generation  depends  on  temperature,  strain  and  strain  rate  [19]. 
To  overcome  this,  a  temperature  dependent  heat  generation  is  derived  as  a  surface  flux  and  given  by,  qtotal  =  corr(T). 
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Here  T  =  T(X,  Y&Z)  and  treated  as  non-uniform  at  the  tool/matrix  interface.  This  model  is  compared  with  the  previous 
linearly  dependent  heat  generation  model  in  which  heat  generation  is  not  dependent  on  temperature  (Figure  6)  [19]. 


xio6 


Figure  6:  Linear  Heat  Generation  Vs  Temperature 
Dependent  Heat  Generation  Model  [19] 

It  was  observed  in  the  probe  tip  that  the  generated  heat  is  same  in  both  the  cases  as  there  is  no  change  in  the 
temperature.  But  in  the  leading  and  trailing  shoulder  region,  q(T)  is  exponentially  increasing.  This  reveals  that  there  is  an 
increase  in  yield  stress  due  to  the  decrease  in  temperature  for  larger  radii. 

CONCLUSIONS 

The  amount  of  heat  generation  in  FSW  depends  on  so  many  parameters  like,  tool  rotation  and  traverse  velocity, 
tool,  tilt  angle,  axial  force,  friction  coefficient,  shear  stress,  duration  of  welding  etc.  The  mutual  dependence  and  influence 
of  all  these  parameters  make  the  analytical  modelling  complicated.  It  is  quite  impossible  to  consider  all  and  normally  the 
most  important  dependences  are  only  considered  in  the  analytical  modelling.  In  order  to  improve  the  accuracy  of  the 
analytical  modelling  it  should  not  rely  only  on  analytical  algebra,  but  also  on  numerical  calculations  and  implementation  of 
experimentally  obtained  variables.  For  example,  implementation  of  a  numerical  flow  model  and  usage  of  experimentally 
estimated  friction  coefficient  and  thermal  conductivity  in  the  analytical  model  gives  a  precise  result. 

Most  of  the  analytical  modelling  considers  a  uniform  heat  generation  in  the  tool  shoulder  metal  interface. 
But  the  frictional  heat  generation  depends  on  the  coefficient  of  friction  between  the  surfaces  in  contact  and  surface 
velocity,  which  are  not  constant  and  but  increasing  gradually  from  the  centre  point  to  the  tool  shoulder  edge  and  this  leads 
to  the  non-uniform  heat  generation  from  the  centre  to  the  shoulder  edge. 

The  difference  in  the  recorded  hardness  in  the  heat  affected  zone  on  the  advancing  side  and  retreating  sides 
reveals  asymmetry  in  the  distribution  of  heat  due  to  the  change  of  direction  of  tool  rotation.  In  order  to  simplify  problems, 
many  thermal  models  assume  the  symmetric  condition  on  both  the  sides  causing  deviation  from  the  experimentally 
measured  values.  In  spite  of  all  these  assumptions  and  its  drawbacks  in  the  analytical  modelling,  it  helps  to  have  a  clear 
physical  understanding  of  the  process  and  to  predict  the  temperature  profile  on  various  stages  so  that  the  process 
parameters  can  be  easily  optimised. 
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